ChinaXiv 合 作 期 刊 


第 40 卷 ”第 5 期 干 早 区 人 研究 Vol.40 No.5 
2023 年 5 月 ARID ZONE RESEARCH May 2023 


长 江 源 区 1980 一 2020 年 水 沙 变化 规律 
Wd o, Aga, a xs 


(1. 中 国 科 学 院 教育 部 水 土 保持 与 生态 环境 研究 中 心 ,陕西 杨凌 712100; 2. 中 国 科 学 院 水 利 部 水 土 保 持 研究 所 ， 
陕西 杨凌 712100; 3. 中 国 科 学 院 大 学 ,北京 ”100049; 4. 长江 水利 委员 会 长 江 科 学 院 , 湖 北 武汉 430010) 


摘 要: 基于 长 江 源 区 1980 一 2020 年 12 个 站 点 日 降水 及 直 门 达 站 年 径流 及 泥 沙 等 数据 ,利用 中 国土 壤 流 失 方程 、 
俩 最 小 二 乘 结构 方程 和 累积 距 平 等 方法 ,分 析 源 区 土壤 侵蚀 .河道 水 沙 变化 及 长 时 间 尺 度 下 与 降水 的 耦合 作用 机 
制 。 结 果 表 明 :(1) 长 江 源 区 年 土壤 侵蚀 模 数 呈 显著 增加 趋势 (P<0.05) , 且 年 均 土 壤 侵蚀 模 数 为 4.71 t hm? a7 , 18 
烈 及 以 上 强度 侵蚀 主要 分 布 在 源 区 东南 部 。(2) 源 区 径流 量 显著 增加 (P<0.05) , 且 2004 年 为 突变 年 , 含 沙 量 和 输 沙 
量 无 明显 变化 趋势 。(3) 降水 对 径流 和 土壤 侵蚀 均 有 显著 正 回 效 应 ,解释 度 分 别 为 70% .52.9% 。 径 流 和 土壤 侵蚀 对 
输 沙 量变 化 均 具 有 直接 正 向 影响 ,降水 、 径 流 和 土壤 侵蚀 对 输 沙 量变 化 的 综合 解释 度 达 72.5% ,结果 可 为 流域 生态 
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工程 实施 效果 的 研究 与 评价 提供 科学 依据 。 


关键 词 : 中 国土 壤 流失 方程 ， 结构 方 程 模型 ， 气 候 变化 ; 土壤 侵 人 


近年 来 , 随 着 气候 变化 和 人 类 活动 的 加 剧 , 极 
端 气候 事件 频 发 ,使 得 流域 径流 `. 泥 沙 等 水 文 要 素 
及 其 输 移 过 程 发 生 改变 "。 国 内 外 许多 学 者 对 尼 罗 
河 “、 密 西西 比 河 ”、 黄 河 “ 和 长 江 “ 等 重要 河流 开 
展 了 河流 水 沙 情势 变化 规律 及 归 因 分 析 研 究 。 其 
中 ,气候 要 素 中 降水 对 河流 水 沙 的 影响 更 为 显著 ， 
人 类 活动 则 是 通过 水 土 保 持 工 程 等 的 实施 减少 流 
域 水 土 流失 ,进而 影响 水 沙 变 化 “”。 不 可 忽略 的 
是 ,流域 水 土 保持 工程 和 水 库 建设 是 近年 来 造成 河 
流 输 沙 量 下 降 的 重要 原因 ,因此 ,开展 流域 水 土 流 
失 与 水 沙 规律 的 综合 研究 对 流域 水 土 保持 和 泥 沙 
治理 工作 十 分 重要 。 

青藏 高 原 是 亚洲 众多 河流 的 发 源 地 ,被 称 为 “ 亚 
洲 水 塔 王 ,对 全 球 气候 变化 的 响应 极为 敏感 ”。 长 
江 源 区 地 处 青藏 高 原 腹 地 ,是 我 国 重要 的 水 源 涵养 
地 与 生态 安全 屏障 区 ” ,区 内 冻 土 、 沼 泽 、 湖 泊 和 高 
寒 湿 地 等 广泛 分 布 , 生 态 环境 极为 脆弱 ”。 近 年 来 ， 


土地 利用 / 履 盖 和 冰雪 融 水 对 径流 变化 也 有 一 定 调 
节 作用 "I。 水 土 流 失 和 入 河 泥 沙 量 增长 的 主要 原 
因 是 降水 的 增加 和 频繁 的 人 类 活动 中 ,而 长 时 间 序 
列 内 植被 变化 则 与 流域 内 泥 沙 变化 相关 性 不 大 。 
Teng 等 “ 泪 和 Wang 等 “基于 模型 模拟 结果 发 现 降水 
对 源 区 水 土 流失 影响 较 大 , 旦 未 来 水 土 流失 量 将 会 
继续 增加 。 由 此 可 见 ,降水 是 引发 源 区 土壤 侵蚀 和 
水 沙 变化 的 主要 驱动 力 ” ,一 方面 直接 影响 坡 面 土壤 
的 侵蚀 ,进而 影响 泥 沙 过 程 “ 7 , 另 一 方面 通过 影响 
径流 变化 间接 影响 河流 携 沙 能 力 ”…。 尽 管 对 长 江 源 
区 径流 和 泥 沙 变化 及 影响 因素 进行 了 多 项 研究 …"， 
但 是 大 /中 流域 尺度 中 侵蚀 产 沙 与 输 移 的 综合 研究 
仍 较 少 ,同时 ,气候 变化 是 影响 源 区 土壤 侵蚀 与 河 
流水 沙 变化 的 重要 因素 。 因 此 ,有 必要 开展 长 江 源 
区 侵蚀 与 水 沙 变 化 及 其 与 气候 变化 耦合 的 响应 关 
系 研 究 ,以 准确 把 握 长 江上 游 泥 沙 来 源 及 流域 侵蚀 
产 沙 对 降雨 在 长 时 间 尺 度 下 的 响应 关系 ,为 源 区 生 


随 着 全 球 气候 变 暖 的 持续 作用 , 源 区 内 水 土 流失 量 
明显 增加 ,径流 量 和 泥 沙 等 要 素 随 之 发 生变 化 "” , 进 
而 对 下 游 河 道 和 流域 管理 产生 影响 ”。 已 有 研究 表 
HH , 源 区 内 径流 的 季节 性 变化 主要 受降 水 的 影响 ”， 
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态 工程 实施 效果 的 研究 与 评价 提供 科学 依据 。 

此 ,本 文 基 于 长 江 源 区 1980 一 2020 年 降水 、 
径流 和 泥 沙 等 资料 ,研究 长 江 源 区 年 土壤 侵蚀 量 、 
径流 量 和 泥 沙 变 化 特征 ,阐明 长 江 源 区 土壤 侵蚀 变 
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化 规律 ,河道 水 沙 变 化 对 地 表 侵 蚀 的 响应 ,以 及 长 
时 间 尺 度 下 泥 沙 变化 对 侵蚀 和 降水 的 响应 关系 ,以 
期 为 推动 小 流域 科学 治理 和 长 江 流域 高 质量 发 展 
提供 理论 依据 。 


1 数据 与 方法 


1.1 研究 区 概况 

长 江 源 区 为 直 门 达 水 文 站 控制 以 上 的 流域 (图 
1) ,地 理 位置 位 于 90"~98"E 31 ~36"N ,总 面积 约 为 
15.7x10' km。 源 区 地 势 西高 东 低 ,平均 海拔 超过 
4000 m, 太 阳 辐 射 较 强 ,日 温差 大 但 年 温差 小 ” ,多 
年 降水 量 显著 增加 (P<0.05) ,年 均 约 为 370 mm ,水 
系 包 括 通天 河水 系 、 沦 沦 河 水 系 , 呈 慎 状 分 布 。 主 
要 土壤 类 型 为 高 山 草原 土 和 高 山 草 馈 土 ,其 中 高 山 
草 甸 土 分 布 面积 最 大 , 约 占 源 区 总 面积 63% ,植被 
类 型 从 东南 到 西北 依次 分 布 有 高 寒 草 多 V qe AE d 
原 ARMIR , 冻 土 面积 约 占 青 藏 高 原 多 年 冻 土 面 
积 的 20%”。 人 口 以 藏族 居民 为 主 ,人 口 密度 不 到 
1 人 km”, 牧 业 是 区 域 经 济 发 展 的 主要 生产 方式 。 
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图 1 长 江 源 区 示意 图 


Fig. 1 Basic information of the headwaters of the Yangtze River 


1.2 数据 收集 与 来 源 

长 江 源 区 直 门 达 水 文 站 1980 一 2020 年 逐年 径 
流 及 泥 沙 数据 由 青海 省 水 文 局 提供 。 从 国家 气象 
科学 数据 中 心 (http:Wdata.cma.cn) 下载 长 江 源 区 及 
其 周围 共 12 个 气象 站 1980—2020 年 逐日 降水 资 
料 。 长 江 源 区 高 程 数据 (空间 分 辩 率 为 1 km) KI 
源 区 归 一 化 植被 指数 (Normalized Difference Vegeta- 
tion Index,NDVI) 数 据 (1998 一 2018 年 ,空间 分 辩 率 为 
1 km) fll 1980 4E. , 1990 ^E. 1995 4E. . 2000 4 , 2005 4F. , 
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2010 年 .2015 年 和 2020 年 8 期 土地 利用 数据 (图 2) 来 
源 于 资源 环境 科学 数据 中 心 (https://www.resdc.cn)。 
1.3 研究 方法 
1.3.1 土壤 侵蚀 量 利用 中 国土 壤 流 失 方 程 
(CSLE) UH KM 1980 一 2020 年 的 土壤 侵蚀 
Æo TE ArcGIS 10.2 的 文 持 下 ,对 各 个 因子 赋值 著 加 
计算 ,得 到 土壤 侵蚀 量 。 模 型 如 公式 (1) 所 示 : 
A-RXKXLXSXBXExT (1) 

式 中 :4 是 土壤 侵蚀 模 数 (t.hm*…a') ;RR 是 降雨 侵蚀 
JJ IF CMJ* mm hm …h'…a');K 是 土壤 可 蚀 性 (t. 
hn?*h*hm?* MJ! mm ) ;是 坡 长 因子 ,S 是 坡度 因 
子 ,二 者 统称 为 地 形 因 子 ,无 量 纲 ;8 是 生物 措施 因 
子 ;五 是 工程 措施 因子 ;7 是 耕作 措施 因子 ,BE 、7 因 
子 均 无 量 纲 ( 取 值 范 围 为 0~1)。 

C1) 降雨 侵 刨 力 因子 (R) 

基于 章 文 波 等 "建立 的 降雨 侵蚀 力 算 法 模型 
计算 年 降雨 侵蚀 力 因子 值 (图 3)， 


R- ay (2) 


式 中 :R 为 半月 周期 的 降雨 侵蚀 力 (MJ.mm:hm*…h… 
a ');E 为 一 个 半月 周期 (每 个 月 的 前 15 d 为 一 个 半 
月 周期 ,该 月 的 其 余天 数 为 男 一 个 半月 周期 );P 是 
一 个 半月 周期 内 超过 12 mm 的 日 降水 量 (mm) ;a 和 
B 分 别 为 模型 的 待定 参数 ,计算 方程 式 如 公式 (3) 和 
公式 (4): 


a=21.586 x p ^! (3) 
[8- 0.8363 + 8.177/Pd,, + 24.455IPy,, (4) 
式 中 :Pd 代表 日 雨量 超过 12 mm 的 日 降雨 量 
(mm) ; P 和 代表 一 年 内 超过 12 mm 的 日 雨量 的 和 
(mm). 
(2) 土壤 可 刨 性 因子 (K) 
土壤 可 刨 性 因子 (K) 是 评价 土壤 被 降雨 分 离 、 
冲 蚀 和 搬运 难 易 程 度 的 一 项 指标 ”。 源 区 土壤 可 
蚀 性 因子 栅 格 数据 集 是 基于 1:1 000 000 中 国土 壤 
数据 库 .传统 的 “土壤 发 生 分 类 "系统 和 Nomo 图 法 
建立 的 K 值 与 土壤 性 质 之 间 的 Nomo 模 型 得 出 。 长 
江 源 区 1000 m 分 辩 率 土壤 可 蚀 性 天 值 来 源 于 国家 
地 球 系统 科学 数据 共享 服务 平台 (http:/www.geoda- 
ta.cn) 。 利 用 ArcGIS 10.2 进行 重 采 样 和 插值 得 到 本 
研究 土壤 可 蚀 性 因子 分 布 (图 4a)。 
(3) 地 形 因子 (LS) 
基于 DEM 数据 计算 地 形 因子 呈 ,得 到 7 值 分 
布 图 (图 4b)。 方 程 如 下 : 
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图 2 长 江 源 区 1980 一 2020 年 土地 利用 变化 
Fig. 2 Landuse changes in the headwaters of the Yangtze River during 1980-2020 
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图 3 长 江 源 区 1980 一 2020 年 年 降水 量 (aj 和 降雨 侵蚀 力 (b) 变 化 


Fig.3 Variations of annual average precipitation (a) and rainfall erosivity factor (b) in the headwaters 


of Yangtze River during 1980-2020 


LS - (GL "(65.41 sin? 8 +4.56 sin B+ 0.065) (5) 


22:13 


式 中 :7 是 地 形 因 子 ;/ 是 累计 坡 长 (m);B8 是 坡度 
(”);m 是 坡度 变量 , 当 坡 度 < 1 时 ,m=0.2; 当 坡 度 为 
1^-3^ Hf ,m20.3; 当 坡 度 为 3.5 ~4.5 时 ,m=0.4; 当 坡 


度 为 25 时 ,m=0.5。 


蚀 的 


HEZ 
E 


(4) 生物 措施 因子 (8) 
生物 措施 因子 (B) 主 要 反映 植被 覆盖 对 土壤 侵 
响 , 其 赋值 范围 为 0~1, 越 接近 0, 表 示 防 治 越 


好 ,几乎 不 发 生 侵蚀 ; 越 接近 1, 表示 未 采取 任何 措 
施 。 具 体 包 括 退 耕 还 林 、 还 草 等 。 本 研究 基于 8 期 
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图 4 长 江 源 区 K 值 (a) 及 LS 值 (b) 空 间 分 布 
Fig. 4 Spatial distribution of K (a) and LS value (b) 


ERGE , 26 Liu 6 Am Fu 等 2 提出 的 赋值 
原则 ( 表 1) ,利用 ArcGIS 10.2 重 分 类 得 到 B 值 的 空 
间 分 布 (图 5)。 


表 1 不 同 土地 利用 类 型 的 8B 值 
Tab.1 B value of different landuse types 


土地 利用 类 型 B 值 土地 利用 类 型 B 值 
旱地 1 城镇 用 地 0.01 
有 林地 0.01 农村 居民 点 0.025 
灌木 林地 0.1 建设 用 地 0.1 
Vi POR 0.1 水 体 0 
其 他 林地 0.2 其 他 未 利用 地 0 
高 覆盖 草地 0.12 沼泽 地 0.05 
中 覆盖 草地 0.18 

低 覆 盖 草 地 0.32 


(5) 工程 措施 因子 (EZ) 和 耕作 措施 因子 (7) 

工程 措施 因子 (E) 反 映 了 实施 的 工程 措施 对 于 
流域 内 水 土 流 失 的 影响 。 长 江 源 区 居民 将 旱地 改 
为 石 坎 水 平 梯田 ,通过 查 表 可 得 石 坎 水 平 梯田 E 值 
为 0.121。 耕 作 措 施 ( 五 ) 通 过 提高 土壤 抗 蚀 性 ,进而 
防治 土壤 侵蚀 。 长 江 源 区 属于 北部 中 高 原 半 干旱 
喜 凉 作物 一 熟 区 ” ,根据 土地 利用 类 型 图 将 研究 区 
耕地 赋值 为 0.121 和 0.488。 其 他 区 域 分 配 为 1, 得 
到 工程 措施 因子 和 耕作 措施 因子 分 布 。 

由 于 长 江 源 区 土地 利用 变化 较 小 3, 因 此 ,各 
期 土地 利用 可 代表 相 邻 5~10 a 变 化 ,基于 CSLE 模 
型 ,估算 得 到 长 江 源 区 1980 一 2020 年 土壤 侵蚀 模 
数 ,根据 4 土壤 侵蚀 分 类 分 级 标准 》(SL190-2007) 将 
长 江 源 区 土壤 侵蚀 模 数 划分 为 6 个 等 级 (0~$ thm ?* 
a 为 微 度 侵蚀 ,5~25 t hm?-a HPE ph ,25-50 t- 
hm”. a" H P EE f fli, 50-80 t hm? a HRZ, 


80-150 t: hm”. a 为 极 强 烈 侵 蚀 ,>150 thm ^a 为 
剧烈 侵蚀 ) ,由 于 微 度 侵蚀 被 定义 为 小 于 土壤 容许 
流失 量 , 因 此 本 文 仅 统计 轻 度 及 以 上 侵蚀 强度 等 级 
的 面积 占 比 1。 

1.3.2 偏 最 小 二 乘 结构 方程 (PLS-SEM) ” 偏 最 小 二 
乘 结构 方程 (PLS-SEM) 是 一 种 结合 因子 分 析 和 回归 
的 非 参 数 方法 ,通过 迭代 估计 最 大 限度 地 解释 独 
立 潜 变 量 的 变化 。 本 研究 利用 PLS-SEM 分析, 通过 
与 多 个 潜在 变量 的 相关 性 来 解释 降水 .土壤 侵蚀 与 
水 沙 变 化 厢 合 模式 。 本 文 所 提出 的 结构 模型 由 降 
KEP) 径流 量 (O) 、 输 沙 量 (S) 和 土壤 侵蚀 量 (M) 
组 成 ,使 用 SmartPLS 3.0 软 件 绘制 4 个 潜 变 量 的 影响 
路 径 , 通 过 300 次 迭代 执行 PLS 算 法 以 获得 路 径 系 
数 (B), 并 使 用 5000 次 抽样 进行 bootstrapping 以 测试 
路 径 系数 的 显著 性 。 当 模型 拟 合 中 NFISO, 
SRMR<0.12 时 , 表明 PLS-SEM 模型 具有 预测 有 效 
TE ,能 较 强 地 解释 因子 间 相 互 影 响 的 变化 。 潜 变量 
间 路 径 是 否 通 过 检验 的 主要 依据 为 显著 性 P 值 , 当 
P<0.05 时 路 径 检 验 通 过 。 

1.3.3 突变 分 析 方 法 利用 累积 距 平 法 和 滑动 : 检 
验 法 2 判断 长 江 源 区 1980 一 2020 年 径流 量 、 含 沙 量 
和 输 沙 量 序列 的 突变 情况 。 


2 结果 与 分 析 
2.1 土壤 侵蚀 的 变化 特征 


图 6 展示 了 长 江 源 区 1980 一 2020 年 各 侵蚀 强 
度 面 积 占 比 变化 。 如 图 所 示 , 轻 度 侵蚀 面积 约 占 总 
侵蚀 面积 的 70.27% ,强烈 及 以 上 强度 侵蚀 面积 约 占 
总 侵蚀 面积 的 10.27%。 长 江 源 区 1980 一 2020 年 土 
壤 侵 蚀 模 数 呈 显著 增加 趋势 (P<0.05)( 图 7), 多 年 平 
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图 5 不 同年 份 生物 措施 因子 值 空 间 分 布 
Fig. 5 Spatial distribution of biological practice factor at the different years 
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图 6 长 江 源 区 1980 一 2020 年 土壤 侵蚀 强度 面积 占 比 变化 
Fig. 6 Changes of soil loss intensity area proportion in the headwaters of Yangtze River during 1980-2020 


15 ERZO 471 t hm7-a^,. Bk (598.16 t 
hm-2.a ,出 现在 1985 年 ;最 小 值 为 2.13 thm ?- a! , 
现在 2015 年。 土壤 侵蚀 模 数 变化 规律 与 降雨 侵蚀 
力 变 化 规律 一 致 (图 3b)。 

图 8 展示 了 长 江 源 区 1980 年 .1990 年 .1995 年 、 
2000 年 、2005 年 .2010 年 .2015 年 和 2020 年 的 土壤 
侵蚀 强度 的 空间 分 布 状况 。 从 中 发 现 源 区 1980 一 
2020 年 土壤 侵蚀 强度 整体 较 低 ,以 轻 度 侵蚀 为 主 ， 
侵蚀 强度 较 大 的 地 方 主 要 分 布 在 东南 部 ,如 玉树 
县 。 东 南部 降水 量 相 对 较 多 , 且 地 形 多 以 高 山 为 


主 , 坡 面 处 极 容 易 发 生 水 力 侵蚀 ,导致 水 土 流 失 情 
况 相 对 其 他 区 域 更 为 严重 。 如 图 7 和 图 8 所 示 , 长 
江 源 区 土壤 侵蚀 时 空 变化 较为 复杂 。1980 年 源 区 
东南 部 土壤 侵蚀 面积 较 大 且 侵 蚀 强 度 较 高 ,1990 一 
1995 年 侵蚀 强度 与 侵蚀 面积 呈现 减 小 趋势 。 土 壤 
侵蚀 模 数 降低 至 3.82 thm-2.a ,强烈 及 以 上 侵蚀 所 
占 面积 明显 减少 ,土壤 侵蚀 得 到 一 定 遇 制 。 然 而 ， 
2000 年 后 侵蚀 面积 和 侵蚀 强度 均 较 高 且 土 壤 侵蚀 
模 数 逐渐 增加 , 轻 度 侵蚀 面积 下 降 , 中 度 及 以 上 侵 
蚀 面 积 占 比 逐渐 增加 。2010 年 后 侵蚀 量 增加 速度 
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图 7 长 江 源 区 1980 一 2020 年 土壤 侵蚀 模 数 变 化 
Fig.7 Change of annual soil loss modulus in the headwaters 
of Yangtze River during 1980-2020 


有 所 放 组 ,其 中 ,2015 年 土壤 侵蚀 面积 最 小 且 整 体 
侵蚀 强度 最 低 ,这 与 2015 年 降雨 侵蚀 力 较 小 有 关 ， 
值 仅 为 125.5 MJ mm*hm?-*h^-a" CE] 3b). 
2.0 年 径流 量 、 含 沙 量 和 输 沙 量 的 变化 规律 

直 门 达 水 文 站 1980 一 2020 年 年 径流 量 .年 含 沙 
量 和 年 输 沙 量 的 变化 趋势 如 图 9 所 示 , 直 门 达 水 文 
站 多 年 径流 量 显著 增加 (P<0.05) ,多 年 平均 径流 量 
为 142.03x10* m’; 年 径流 量 最 大 值 在 2009 年 ,为 
245.14x10* ,最 小 值 为 81.95x10 m ,出 现在 1997 
年 。 相 比 于 90 年 代 ,21 世 纪 以 来 径流 量 明 显 增加 ， 
增幅 约 为 79.7%。 多 年 含 沙 量 和 输 沙 量 均 无 显著 变 
化 趋势 ,多 年 平均 含 沙 量 为 0.69 kgm’, EK [EE HER 


土壤 侵蚀 模 数 /(t:hm ?a1) 
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9 ， BINH 在 1989 年 ,为 1.41 kg.m-2 ,最 小 值 出 现在 1994 年 ,为 
8 0.22 kg'm-?。 多 年 平均 输 沙 量 为 1063.41x10't6 输 沙 
— 10 a2 Ea yh E TES E 
7 E s 量 与 含 沙 量变 化 较为 一 致 ,最 大 值 均 出 现在 1989 4. 
^ 为 2980x10't, 最 小 值 均 出 现在 1994 年 ,为 203x10+t。 


2.3 kk, HE BPURUZK PP BR A AR RR 

基于 1980 一 2020 年 年 降水 量 (P) .年 土壤 侵蚀 
E CM) .年 径流 量 (0O) 和 年 输 沙 量 ($) 数 据 ,通过 
PLS-SEM 分 析 ,得 到 长 江 源 区 1980 一 2020 年 降水 、 
土壤 侵蚀 与 水 沙 变化 耦合 关系 图 (图 10)。 计 算得 
知 , i708 NFI-0.982 (NFI>0.9) , SRMR-0.03 (SRMR« 
0.12) ,表明 模型 对 潜 变 量 的 拟 合 较 好 。1980 一 2020 
年 长 江 源 区 各 潜 变 量 之 间 的 主要 关系 如 下 :(1) 以 
降水 为 驱动 因素 ,降水 对 土壤 侵 刨 和 河道 径流 均 具 
有 直接 的 显著 正和 加 效应 [| 6=0.727(P<0.01) 和 B= 
0.837 (P«0.01) ] ,并 进一步 通过 径流 和 侵蚀 产 沙 间 
接 影响 输 沙 量变 化 ,(2) 径流 和 土壤 侵蚀 对 输 沙 量 
变化 具有 直接 正 问 效应 [| 6=0.885 (P<0.01) 和 B= 
0.062 ] ,径流 对 输 沙 量变 化 的 贡献 度 明 显 大 于 土壤 
侵蚀 。 通 过 间接 效应 分 析 发 现 , 降 水 对 输 沙 量具 有 
显著 正身 影响 。 降 水、 径流 和 土壤 侵蚀 对 输 沙 量变 
化 的 解释 度 达 72.5% , 其 影响 路 径 主要 为 降水 一 径 
流 一 输 沙 ( 表 2)。 


3 讨论 


3.1 长 江 源 区 土壤 侵蚀 对 泥 沙 变化 的 影响 
本 研究 利用 CSLE 模 型 估算 1980 一 2020 年 长 江 
源 区 多 年 平均 土壤 侵蚀 模 数 的 结果 与 Liu E IRR 


1980 年 Hu 19904 < 1995 年 2000 年 À 
E 2005 年 | 2010 年 | 2015 年 2020 年 


0 125km 


图 8 长 江 源 区 1980 一 2020 年 土壤 侵蚀 模 数 空间 分 布 
Fig. 8 Spatial distribution of soil loss modulus in the headwaters of Yangtze River during 1980-2020 


202305.00252v1 


chinaXiv 


rk, Tou 


径流 量 /108 m 


Arb Rem) 


y=4.07x—7092.68 
R?-0.006 


输 沙 量 /104t 


1980 1985 1990 1995 2000 2005 2010 2015 2020 
年 份 


图 9 长 江 源 区 1980 一 2020 年 径流 量 (a)、 含 沙 量 (b) 和 输 沙 
量 (c) 的 变化 趋势 

Fig. 9 Variations of runoff (a), suspend sediment concentrate 

(b), and sediment load (c) in the headwaters of Yangtze River 
during 1980-2020 


注 :P 为 年 降水 量 ,为 年 土壤 侵蚀 量 ,Q 为 年 径流 量 ,5 为 年 输 沙 
量 ;** 表 示 0.01 水 平 上 相关 ,* 表 示 0.05 水平 上 相关 。 下 同 。 

图 10 长 江 源 区 1980 一 2020 年 降水 .土壤 侵蚀 与 水 沙 变化 

耦合 关系 


Fig. 10 Coupling relationship between precipitation, soil 


loss, runoff, and sediment load in the headwaters of Yangtze 
River during 1980-2020 


梦 美 等 ”的 分 析 结 果 较 为 一 致 。 但 是 ,显著 高 于 
Teng 等 "人 舍 算 的 2.76t.hm*…a'。 产 生 差 异 的 主要 
原因 是 计算 相关 因子 值 的 数据 来 源 及 估算 方法 不 
同 ,特别 是 RR 因子 和 kK 因子。 由 于 所 选 气象 站 点 数 
量 和 位 置 的 不 同 ,计算 得 到 的 降雨 侵蚀 力 具 有 一 定 
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R2 1980 一 2020 年 长 江 源 区 降水 、 土 壤 侵 蚀 、 径 流 
和 输 沙 量 路 径 效 应 
Tab. 2 Route effects of precipitation, soil loss, runoff, and 


sediment load in the headwaters of Yangtze River during 


1980-2020 
路 径 直接 效应 间接 效应 总 效应 
P>M 0.727 - 0.727 
PpP—0 0.837 - 0.837 
PS - 0.696 0.696 
0 一 5 0.885 - 0.885 
M>S 0.062 - 0.062 


UE: P KERKE , M NFER , Q 为 年 径流 量 ,$ 为 年 输 沙 量 。 


差异 ,本 人 研究 中 年 均 R 值 (269.25 MJ .mm hm.h 
a ) 5T xr ? R mE RR fc ph 7J v5. El (25.72~ 
448.35 MJ .mm.hm-2.h 7a) 略 有 不 同 , 但 在 可 接受 
变化 范围 内 。 同 时 ,区 域 土 壤 侵蚀 量 和 降雨 侵蚀 力 
的 空间 变化 均 为 由 西北 向 东南 增加 。 另 外 ,本文 对 
K 值 栅 格 数据 集 进行 重 采样 及 空间 插值 后 ,发现 源 
区 侵蚀 主要 集中 在 草 甸 土 ,与 全 国 水 土 流失 与 生态 
安全 考察 结果 一 致 ,这 也 使 得 降雨 人 渗 量 降低 进而 
导致 径流 量 增加 ” 。 

源 区 多 年 平均 土壤 侵蚀 模 数 呈 显 车 增加 趋势 ， 
然而 ,多 年 输 沙 量 和 含 沙 量 无 明显 变化 趋势 。 一 方 
面 ,可 能 是 源 区 水 力 侵蚀 主要 分 布 在 源 区 东南 部 ， 
该 区 域 以 高 山 为 主 , 且 7ZS 值 较 高 ,进而 使 泥 沙 容易 
在 坡 面 洲 积 ; 另 一 方面 ,虽然 新 的 土壤 侵蚀 仍 在 增 
加 ,但 源 区 水 土 保持 工作 取得 了 一 定 成 效 , 源 区 输 
沙 量 整体 处 于 一 个 相对 平衡 的 状态 ”。 另 外 ,20 世 
纪 80 年 代 源 区 人 为 无 序 开采 沙 金 导致 源 区 的 侵 刨 
面积 和 侵蚀 强度 都 较 大 ” ,整体 侵蚀 情况 较为 严 
重 ;20 世 纪 90 年 代 源 区 植被 生态 发 展 趋势 良好 , 没 
有 出 现 退 化 区 域 ' 引 ,生态 环境 总 体 趋势 良好 ,因此 ， 
水 土 流失 状况 得 到 明显 改善 ;2000 一 2004 年 , 近 
70% 的 地 区 出 现 植被 轻微 退化 的 现象 ”, 且 长 江 源 
区 降雨 侵蚀 力 增加 ,导致 影响 水 土 流失 的 外 部 动力 
条 件 增强 ” ,进一步 加 剧 了 水 土 流失 状况 。 但 是 ， 
得 益 于 近 些 年 三 江 源 自然 生态 保护 规划 的 实施 55 , 
2010 年 后 侵蚀 增加 速度 有 所 放 组 (网 7)。 

尽管 近年 来 源 区 进行 了 大 量 的 生态 保护 与 恢 
复工 作 ,但 是 对 比 前 期 水 土 流失 遥感 调查 数据 发 
现 ,土壤 侵蚀 仍 在 加 剧 ,对 生态 环境 的 威胁 依然 存 
在 ,长 江 源 区 的 水 土 保护 及 生态 修复 工作 仍 需 继 续 
开展 。 
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3.2 不 同时 期 降水 、 土 壤 侵 蚀 与 水 沙 的 响应 关系 的 
差异 

由 1980—2020 年 长 江 源 区 年 均 径 流量 累积 距 
平 图 可 知 ( 图 11),2004 年 是 年 径流 量 累 积 距 平 值 发 
生 显 著 变化 的 一 年 ,在 2004 年 达到 最 低 ,随后 升 
高 。 结 合 滑动 1 检验 法 ,进一步 确定 了 2004 年 为 径 
流量 的 突变 点 。 年 输 沙 量 累 积 距 平 值 在 1988 年 和 
2006 年 发 生变 化 ,但 两 个 突变 点 均 未 通过 检验 ,年 
含 沙 量 在 1994 年 发 生 的 突变 也 未 通过 检验 。 因 此 ， 
根据 径流 量 突变 点 将 水 文 序列 划分 为 1980—2004 
年 .2005 一 2020 年 两 个 阶段 。 


—e— 1980—20044£ 
——- 2005—20204E 


径流 量 累积 距 平 /10 m? 


1980 1985 1990 1995 2000 2005 2010 2015 2020 
年 份 


图 11 长 江 源 区 1980 一 2020 年 径流 量 累 积 距 平 变化 


Fig. 11 Cumulative anomaly changes ofrunoff in the 
headwaters of Yangtze River during 1980-2020 


源 区 直 门 达 水 文 站 输 沙 量 较为 稳定 ,径流 量 在 
2004 年 出 现 抛 点 ,与 刘 彦 等 ”和 张 永 勇 等 ”的 研究 
结果 一 致 。 径 流量 发 生 突变 的 原因 主要 是 21 世纪 
初 降水 的 增加 与 三 江 源 生态 保护 项 目的 实施 ”。 
此 ,将 2004 年 作为 分 界 点 ,对 1980 一 2004 年 、 


(a) 1980—20044£ 


0.690 


P 0.981” 
河道 径流 T 
土壤 侵蚀 


降水 c0 — 输 沙 


2005 一 2020 年 两 个 时 期 进行 耦合 关系 分 析 可 知 (图 
12) ,降水 对 径流 量变 化 的 解释 度 上 升 [| B=0.690(P< 
0.01) 和 B=0.869(P<0.01) ] ,结合 图 10 分析,3 个 时 期 
中 模型 均 显 示 降 水 与 河道 径流 呈 显 著 正 相关 , 且 降 
水 对 径流 变化 贡献 度 达 70% 左 右 ,降水 是 径流 量变 
化 的 主导 因素 ,与 李 凯 等 ”的 研究 结果 较为 一 致 。 

通过 对 图 10 和 图 12 中 3 个 时 期 耦合 关系 结 
对 比 发 现 降水 对 土壤 侵蚀 的 正 向 影响 也 略 有 上 升 ， 
且 均 达到 显著 水 平 。 结 合 降 雨 侵蚀 力 ( 图 3b) 与 十 
壤 侵 蚀 变 化 ,认为 降水 是 引起 源 区 侵蚀 的 主要 因 
素 。 但 是 , 源 区 仍 面临 由 于 人 类 活动 导致 草场 退化 
及 土地 沙化 的 严峻 形势 ” ,未 来 有 可 能 会 加 剧 源 区 
侵蚀 并 影响 源 区 径流 与 泥 沙 变化 。 

目前 大 多 数 研 究 认 为 区 域 土壤 侵蚀 对 输 沙 量 
有 较 大 影响 ” ,然而 ,本 文 研究 结果 显示 土壤 侵 伺 
与 输 沙 量变 化 没有 显著 相关 关系 , 源 区 土壤 侵蚀 与 
含 沙 量 、 输 沙 量变 化 趋势 也 明显 不 同 。 一 方面 ,水 
文 特性 和 流域 的 形态 等 因素 也 会 引起 输 沙 量变 化 ; 
男 一 方面 由 于 源 区 侵蚀 泥 沙 在 坡 面 淤积 且 源 区 内 
河 湖 、 湿 地 等 存储 效应 的 影响 ,从 而 导致 输 沙 量 
变化 较为 稳定 ,所 以 源 区 输 沙 量 的 变化 与 侵蚀 并 没 
有 表现 出 必然 联系 “|。 
3.3 土地 利用 及 植被 覆盖 变化 对 流域 水 沙 的 影响 

土地 利用 / 窗 被 变化 对 流域 水 文 活动 起 着 重要 
的 调节 作用 ”, 其 对 流域 产 流产 沙 的 影响 是 流域 水 
文学 和 生态 水 文学 研究 的 关键 问题 之 一 。 

21 世纪 初期 长 江 源 区 人 类 扰动 较 显 闭 ,2010 年 
后 较为 平稳 后。 在 本 研究 时 段 内 (1980 一 2020 年 ) 
源 区 草地 多 年 平均 占 比 约 在 68%~70% ,其 次 是 未 利 
用 地 和 水 域 , 近 40 a 土地 利用 格局 相对 稳定 (图 
2)。 基 于 1998 一 2009 年 .2010 一 2018 年 土地 利用 和 
1 km 空间 分 辩 率 的 归 一 化 植被 指数 (NDVIT) 数 据 与 


(b) 2005—20204Ẹ 
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输 沙 
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图 12 长 江 源 区 1980—2004 EFI 2005—2020 EKEk .土壤 侵蚀 与 水 沙 变化 耦合 关系 
Fig. 12 Coupling relationship between precipitation, soil loss, runoff, and sediment load in the headwaters of Yangtze River during 
1980-2004 and 2005-2020 
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水 沙 进行 相关 性 分 析 发 现 ( 表 3) ,土地 利用 和 植被 
指数 变化 与 水 沙 变 化 整体 几乎 呈现 不 显著 相关 关 
系 ,与 郭 帅 等 2 和 Zhang 等 5 的 研究 结果 相同 。 同 
时 ,Naveed 等 ”和 于 等 ”认为 长 江 源 区 径流 量变 化 
主要 由 气候 变化 引起 ,与 土地 利用 相 比 ,气候 变化 
的 影响 要 大 得 多 。 男 外 ,在 本 研究 中 土地 利用 的 变 
化 已 吉 括 于 土壤 侵蚀 的 计算 中 ,因此 ,本文 重点 探 
讨 了 降水 对 水 沙 变 化 的 影响 ,不 再 单独 讨论 土地 利 
用 和 植被 覆盖 变化 对 流域 侵蚀 与 水 沙 变化 的 影 
响 。 随 着 人 类 活动 的 增加 ,加 之 气候 变化 的 影响 和 
长 江 源 区 各 项 生态 项 目的 不 断 推进 ,未 来 长 江 源 区 
土地 利用 和 植被 覆盖 因子 的 变化 不 可 忽视 ,应 注重 
未 来 源 区 生态 保护 及 水 源 涵养 |。 


表 3 1998 一 2009 年 .2010 一 2018 年 土地 利用 、 植 被 覆盖 
变化 与 径流 、 泥 沙 相 关 性 分 析 
Tab.3 Correlation analysis of landuse, NDVI, runoff, 


and sediment load from 1998-2009, 2010-2018 


时 期 NDVI ”草地 水 体 ”未 利用 地 

1998 一 2009 径流 量 0.538 -0.146 0.403 -0.028 

i 输 沙 量 ^ 0.426  -0.068 0145 . -0.184 
土壤 侵蚀 -0.300 - - - 

2010—2018 径流 量 ^ 0.614 . 0.302 -0.302 0.302 

a 输 沙 量 ^ 0.619 0.135 -0.135 0.135 
土壤 侵蚀 0.740 - - - 


注 :* 表 示 0.05 水 平 相关 。 


4 结论 


本 文 基 于 中 国土 壤 流失 方程 和 结构 方程 模型 
分 析 长 江 源 区 水 土 流失 与 水 沙 变化 特征 及 其 与 降 
水 的 啊 应 关系 ,得 出 以 下 结论 : 

(1) 长 江 源 区 1980 一 2020 年 土壤 侵蚀 强度 整 
体 较 低 ,以 轻 度 侵蚀 为 主 , 约 占 总 侵蚀 面积 的 
70.27% ,年 均 土壤 侵蚀 模 数 为 4.71 thm”. a”, Œ Œ 
著 增 加 趋势 (P<0.05) ,但 在 2000 年 左右 增加 趋势 最 
为 明显 ,主要 是 由 于 降水 的 影响 。 侵 蚀 强度 分 布 具 
有 明显 空间 异 质 性 ,强度 较 大 的 地 方 主要 分 布 在 源 

(2) 1980 一 2020 年 直 门 达 水 文 站 多 年 径流 量 呈 
显著 增加 趋势 ,年 均 径 流量 为 142.03x10: m? , 输 沙 量 
和 含 沙 量 均 无 明显 变化 趋势 ,水 沙 变 化 具有 较 好 的 
一 致 性 ,径流 量 在 2004 年 发 生 突变 , 突变 后 呈现 增 
加 趋势 。 

(3) 降水 对 土壤 侵蚀 .径流 和 输 沙 均 具 有 显著 
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正 向 效应 , 且 降 水 是 土壤 侵蚀 发 生 的 主导 因素 , 降 
水 主要 通过 影响 径流 对 输 沙 变化 产生 正 向 影响 。 
未 来 应 及 时 关注 降水 带 来 的 土壤 侵蚀 影响 ,加 强 未 
来 源 区 生态 保护 和 恢复 工作 。 
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Abstract: The changes in soil loss, runoff, and sediment load in the headwaters of the Yangtze River and their 
correlation were studied based on the data of daily precipitation at 12 stations, and annual runoff and sediment at 
Zhimenda station in the headwaters of this river from 1980 to 2020. (1) The annual soil erosion modulus in the 
source area showed a significantly increasing trend (P«0.05), and the average annual soil erosion modulus was 
4.71 t-hm^-a''. The erosion intensity and higher grades of erosion were mainly distributed in the southeast of the 
headwaters of the Yangtze River. (2) The runoff increased significantly (P«0.05), however, there was an abrupt 
change in 2004 when the suspended sediment concentrate and sediment load showed no significant change. (3) 
Precipitation had a significant positive effect on runoff and soil loss at 70% and 52.995, respectively. Both runoff 
and soil loss showed direct positive effects on the sediment load. Precipitation, runoff, and soil loss explained 
72.5% of the variation in sediment load. During this period, changes in precipitation had significant effects on soil 
loss and sediment in the headwaters of the Yangtze River. The results can provide a scientific basis for the research 
and evaluation of the effects of implementing ecological engineering in the headwaters of the Yangtze River. 
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